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CRITICAL  SHEAR  STRESS  OF  INFINITELY  LONG-,  SIMPLY 
SUPPORTED  PLATE  WITH  TRANSVERSE  STIFFENERS 
By  Manuel  Stein  and  Robert  W.  Frail  ch 


SUMMARY 


A  theoretical  solution  is  given  for  the  critical  shear  stress  of 
an  infinitely  long,  simply  supported,  flat  plate  with  identical,  equally 
spaced,  transverse  stiffeners  of  zero  torsional  stiffness.  Results  are 
obtained  by  means  of  the  Lagrangian  multiplier  method  and  are  presented 
in  the  form  of  design  charts.  Experimental  results  are  included  and 
are  found  to  be  in  good  agreement  with  the  theoretical  results. 


INTRODUCTION 


The  design  of  shear  web  beams  and  nonwrihkling  skin  surfaces 
requires  a  knowledge  of  the  critical  shear  stress  of  stiffened  plates. 
The  purpose  of  the  present  paper  is  to  give  the  theoretical  critical 
shear  stress  of  an  infinitely  long,  simply  supported,  flat  plate  rein¬ 
forced  with  identical,  equally  spaced,  transverse  stiffeners. 

The  results  are  found  by  means  of  the  Lagrangian  multiplier  method. 
The  stiffeners  are  assumed  to  have  bending  stiffness  but  no  torsional 
stiffness  and  are  assumed  to  be  concentrated  along  transverse  lines  in 
the  middle  plane  of  the  plate.  The  assumption  that  the  stiffeners  have 
no  torsional  stiffness  applies  with  little  error  in  the  case  of  many 
open  section  stiffeners.  The  assumption  that  the  stiffeners  are  con¬ 
centrated  along  transverse  lines  in  the  middle  plane  of  the  plate  is 
applicable  whenever  the  width  of  the  attached  flange  is  small  in  com¬ 
parison  with  the  stiffener  spacing. 

The  theoretical  analysis  of  the  problem  is  given  in  the  appendixes. 
For  completeness,  an  energy  solution  for  the  plate  with  relatively  weak 
stiffeners  is  given  in  appendix  A.  The  solution  for  a  plate  with 
stiffeners  of  intermediate  or  higher  bending  stiffness  is  given  in 
appendix  B.  The  results  are  presented  in  the  form  of  nondimens ional 
curves  which  cover  the  complete  range  of  stiffener  stiffness  and  various 
stiffener  spacings  and  in  a  table  giving  values  from  which  the  curves 
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vere  dravn  (table  i) .  Experimental-  results  are  presented  for  20  panels. 
Comparison  of  these  results  with  the  present  theory  indicates  good 
agreement  between  theory  and  experiment. 


SYMBOLS 


t 

b 

d 

b/d 

D 

E 

P 

E 

I 

P 

El 

Dd 

X 

v 

(Vs)l 


critical  shear  stress 


critical  shea3>-^tress  coefficient 


thickness  of  the  plate 
width  of  plate 
stiffener  spacing 
panel  aspect  ratio 

flexural  stiffness  of  the  plate 


E  t3 
P 


12(1  -  H2) 


Young1  s  modulus  for  plate 


Young1 s  modulus  for  stiffener 
effective  moment  of  inertia  of  stiffener 
Poisson1 s  ratio  for  material 


ratio  of  stiffener  stiffness  to  plate  stiffness 


half  wave  length  of  buckles 
deflection  of  the  plate 
deflection  of  the  ifth  stiffener 

reference  axes 


m,  n, 
r,  i 


integers 
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V 


coefficients  of  deflection  function 


undetermined  Lagrangian  multipliers 

internal  energy  of  tending  of  the  plate 
internal  energy  of  tending  of  stiffeners 

external  work  of  the  stresses 


BACKGROUND 


The  problem  of  the  buckling  of  stiffened  plates  in  shear  has  teen 
treated  by  many  authors  by  the  use  of  both  theoretical  and  semi- 
empirical  methods.  In  1930  Schmieden  (reference  1)  solved  the  differ¬ 
ential  equation  for  an  infinitely  long  plate  stiffened  by  closely 
spaced  transverse  stiffeners  (equivalent  to  orthotropic  plate)  and 
found  exact  stability  criterions  for  shear  buckling  of  plates  with 
simply  supported  edges  and  with  clamped  edges.  By  making  certain 
simplifying  modifications  of  the  stability  criterions*  Schmieden 
obtained  approximate  values  of  the  critical  shear  stresses.  Later  in 
1930  Seydel  (reference  2)  obtained  exact  solutions  for  infinitely  long 
orthotropic  plates  with  simply  supported  or  clamped  edges.  With  the 
use  of  the  proper  parameters  Seydel* s  results  can  be  readily  applied 
to  plate— stiffener  combinations.  The  values  of  the  stresses  obtained 
from  Schmieden1 s  theory  lie  slightly  below  the  exact  values  of  Seydel. 
In  1947  T.  K.  Wang  (reference  3)  used  the  energy  method  to  obtain  an 
approximate  solution  for  plate— stiffener  combinations  with  simply 
supported  edges.  Wangfs  results  lie  above  the  exact  values  of  Seydel. 
All  the  foregoing  solutions  are  applicable  only  to  the  case  of  weak 
stiffeners*  where  the  stiffening  effect  of  the  stiffeners  can  be 
considered  to  be  uniformly  distributed  over  the  plate. 

Solutions  are  also  available  for  plates  reinforced  by  rigid 
stiffeners.  In  193&  Timoshenko  (reference  4)  treated  the  case  of 
simply  supported  rectangular  plates  reinforced  with  one  or  two 
stiffeners.  By  means  of  th9  energy  method  Timoshenko  found  the 
stiffener  flexural  rigidity  necessary  to  prevent  buckling  across 
stiffeners  with  the  conservative  assumption  that  the  stiffeners  act  as 
simple  supports.  In  1948*  Budiansky*  Conner*  and  Stein  (reference  5) 
found  the  critical  shear  stress  for  an  infinitely  long*  clamped  plate 
divided  into  square  panels  by  nondeflecting  intermediate  supports  which 


a  *  b  * 
n*  n* 

^mn 

V 
Vs 


T 
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correspond  to  rigid  stiffeners.  They  also  considered  the  case  of  a 
plate  of  infinite  length  and  width  having  nondeflecting  intermediate 
supports  that  form  an  array  of  square  panels. 

Kuhn  has  written  a  number  of  papers  on  related  subjects  in  which 
he  presents  semiempirical  results  for  the  critical  shear  stress  of 
stiffened  plates.  (See,  for  example,  reference  6.) 

The  available  theoretical  solutions  treat  the  relatively  unim¬ 
portant  case  of  weak  or  closely  spaced  stiffeners  and  the  case  of 
rigid  stiffeners  that  divide  a  plate  into  square  panels.  None  of  the 
theoretical  solutions  presents  results  for  the  practical  range  of 
intermediate  stiffener  stiffness  and  very  little  theory  is  presented 
for  the  practical  range  of  spacing  of  rigid  stiffeners.  Also,  it  is 
felt  that  the  semiempirical  results  for  transverse  stiffened  plates 
cannot  be  extended  to  all. stiffener  spacings  and  stiffnesses  without 
a  sound  theoretical  basis.  The  theoretical  results  of  the  present 
paper  cover  the  complete  range  of  stiffener  stiffness  and  the  practical 
range  of  stiffener  spacing. 


RESULTS  AND  DISCUSSION 


The  critical  shear  stress  for  a  plate— stiffener  combination 
is  given  by  the  formula 


r  =  k 


8  b2t 


Curves  are  presented  in  figure  1  giving  corresponding  values  of  kg 

and  the  stiffness  parameter  ^  for  simply  supported,  transversely 

Da 

stiffened  plates  with  panel  aspect  ratios  of  1,  2,  and  5.  These 
results  are  replotted  in  logarithmic  form  in  figure  2  for  comparison 
with  experimental  results. 

The  points  of  discontinuity  of  the  slopes  in  the  curves  of  figure  1 
represent  changes  in  buckle  patterns.  The  present  results  for  an  ortho¬ 
tropic  plate  agree  with  the  exact  results  of  reference  2.  The  deri¬ 
vation  of  the  buckling  criterion  for  an  orthotropic  plate  (a  plate 
stiffened  by  stiffeners  of  low  bending  stiffness)  is  given  in  appendix  A 
The  derivation  of  the  buckling  criterion  for  plates  stiffened  by 
■stiffeners  of  higher  bending  stiffness  is  given  in  appendix  B. 
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In  previous  solutions*  values  of  kg  -were  found  by  using  the 
orthotropic— plate  curve  and  a  cut-off  at  the  value  of*  kQ  for  simply 

supported  panels.  (See  fig.  1.)  These  figures  show  that  the  present 
solution  yields  values  of  kg  that  are  considerably  below  those  given 

by  the  orthotropic— plate  curve  in  the  intermediate  range  of  stiffener 
stiffness.  Also*  the  present  solution  for  more  rigid  stiffeners  yields 
a  curve  that  is  higher  than  the  cut-off*  which  is  obtained  by  assuming 
the  stiffeners  to  have  the  effect  of  simple  supports.  Since  the  conti¬ 
nuity  of  the  plate  across  the  stiffeners  of  higher  bending  stiffness 
certainly  adds  a  constraint  to  the  plate*  a  higher  buckling  stress  than 
that  corresponding  to  a  simply  supported  edge  is  obtained. 

In  figure  2*  experimental  results  are  compared  with  the  theoretical 
curves.  These  results  are  from  two  sources.  The  first  set  of  experi¬ 
mental  data  is  taken  from  MCA  tests  on  shear  webs  of  24S-T  aluminum 
alloy  attached  to  torsion  boxes.  Drawings  of  a  shear  web  and  torsion 
box  and  the  method  of  loading  are  given  in  reference  7-  Buckling  loads 
were  obtained  from  the  stiffener  load-deflection  curves  which  were 
taken  from  the  original  data.  Each  of  the  buckling  loads  given  in  the 
present  paper  is  the  average  load  at  which  the  stiffeners  start  to 
deflect.  The  properties  of  the  specimens  and  the  buckling  data  are 
given  in  table  II. 

The  second  set  of  experimental  data  is  taken  from  MCA  tests  on 
thick  web  beams  described  in  reference  3.  The  beams  were  made  of 
24S-T  aluminum  alloy  with  heavy  flanges  and  with  joggled  stiffeners 
riveted  to  the  flanges.  The  open  spaces  in  the  joggles  were  filled 
with  soft  metal.  A  picture  of  a  failed  beam  is  shown  in  figure  3. 

The  load  was  applied  at  the  center  and  the  reactions  were  at  the 
ends  of  the  beams.  Lateral  deflections  were  prevented  by  lateral 
supports.  The  load*  when  strain  was  first  observed  in  the  stiffeners* 
was  taken  as  the  buckling  load.  The  properties  of  the  specimens  and 
the  buckling  data  are  given  in  table  III. 

The  stiffener  spacings  for  the  test  results  are  not  the  same  as 
those  for  the  theoretical  results.  All  the  test  results  fall  in  the 
expected  regions  among  the  theoretical  curves.  Only  the  group  of  test 

results  for  which  -  =  2.4  fall  in  the  range  which  serves  to  verify 

d 

the  present  theory  over  previous  theory  which  considered  the  orthotropic— 
plate  curve  to  hold  up  to  the  cut-off  at  which  the  stiffeners  are 
assumed  to  act  as  simple  supports.  The  other  groups  of  test  results 
agree  with  the  present  theory*  but  they  do  not  cover  the  range  in  which 
an  appreciable  difference  exists  between  the  present  theory  and  previous 
theory.  More  experimental  results  are  required  to  confirm  the  present 
theory  fully. 
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CONCLUDING  REMARKS 


Charts  ar©  presented  from  which  th©  th©or©tical  critical  sh©ar 
stresses  can  he  obtained  for  infinitely  long,  simply  supported  plates 
stiffened  with  identical,  equally  spaced,  transverse  stiffeners  of  zero 
torsional  stiffness.  The  theoretical  results  are  "based  on  the  Lagrangian 
multiplier  method.  Previous  theory  considered  the  orthotropic  curve 
to  hold  up  to  a  cut-off  value  corresponding  to  the  stiffener  stiffness 
at  which  the  buckling  load  was  equal  to  the  buckling  load  of  a  simply 
supported  panel  the  size  of  each  bay.  Comparison  of  the  present  theory 
and  previous  theory  shows  that  previous  theory  gives  unconservative 
results  for  stiffeners  of  intermediate  stiffness  and  conservative 
results  for  stiffeners  of  high  stiffness.  Test  results  of  20  panels 
are  presented  which  are  in  good  agreement  with  the  present  theory.  For 
a  conclusive  check  additional  test  results  are  required. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 

Langley  Air  Force  Base,  Va.,  January  28,  1949 
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APEENDIX  A 

THEORETICAL  SOLUTION  OF  CRITICAL  SHEAR  STRESS  OF  ELATES 
WITH  TRANSVERSE  STIFFENERS  CF  LOW  BENDING  STIFFNESS 


If  the  stiffener  handing  stiffness  is  low  and  the  stiffeners  are 
fairly  closely  spaced,  the  buckle  pattern  may  he  considered  independent 
of  the  stiffener  spacing,  and  the  plate  stiffener  combination  can  then 
he  analyzed  as  a  plate  with  different  bending  properties  in  each 
direction,  that  is,  an  orthotropic  plate.  In  this  appendix  buckling  in 
shear  of  an  orthotropic  plate  is  analyzed  by  means  of  the  energy  method. 

The  buckling  configuration  of  the  plate  shewn  in  figure  4  is 
represented  by  the  trigonometric  series 


w 


m 


which  satisfies  the  boundary  conditions  of  simple  support  term  by  term. 

The  internal  bending  energy  of  the  plate  V,  the  internal  bending  energy 

of  the  stiffeners  Va,  and  the  external  work  of  the  shear  stresses  T 

s 

are  given  by  the  expressions 
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Substitution  of  the  expansion  for  w  (equation  (Al))  into  these 
energy  integrals  gives 


n=l*3*  •  •  •  <3.— • 


Then 


n=l,3^-*»  q=2^4j... 
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where 


k 


s 


According  to  the  energy  method  the  potential  energy  (Y  +  V  —  T) 
nrnat  he  minimized  with  respect  to  the  unknown  coefficients  and  b 

By  minimizing  (V  +  Vg  -  T)  with  respect  to  the  coefficients  a  and  b  , 
the  following  set  of  equations  is  obtained:  D 


Sn 


b2  2\2  ^  El 

—  +  nd)  +  n  — 

/  Dd 


8bk£ 

Jtl 


2 


nq 


r  =  0 


1#  3f  •  • 


.  V-"2) 

(n=2,4,6,...) 


(A4) 


n 


r.o 

k  +  „2N)  +  n*  51 
v2  J  Bd 


Sb^ 


2 


q=2,4,... 


a  _ 2S.  ,  0 

5(»2-*2) 

(n=l_,  3*5* • • • ) 


(A5) 


The  coefficients  a^  can  be  found  in  terms  of  br  from 
equation  (A4).  Substitution  of  the  resulting  expression  for  a  in 
equation  (A5)  results  in  the  following  equations: 


b 


n 


+ 


(A6) 


(n— 1*  3 >5^ • • • ) 
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A  solution  to  equations  (A6)  exists  if  the  following  stability  determi¬ 
nant  vanishes: 


Cn  ci3  ci5  •  •  • 

C31  C33  c35 


(where  n  ^  r) 


A  solution  including  all  the  a  *s  and  b  can  he  obtained  by 

n  1 

setting  equal  to  zero  the  first  approximation  of  the  determinant 
equation  (A7) 
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I 

Similarly  the  second  approximation  includes  all  the  a^1 s,  h^, 
0  and  l»2 


C11  c33 


=  0 


Higher  approximations  are  found  in  a  similar  manner.  A  second  approxi¬ 
mation  was  found  to  give  satisfactory  results.  For  a  given  approxi¬ 
mation  it  is  necessary  to  try  values  of  b/X  and  find  the  corre¬ 
sponding  values  of  ka  until  a  minimum  value  of  ks  with  respect 

ET 

to  b/X  is  found  for  each  — .  The  results  are  given  in  table  I  and 

Dd 

in  figure  1# 
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APPENDIX  B 

THEORETICAL  SOLUTION  OF  CRITICAL  SHEAR  STRESS  OE  PLATES 
WITH  TRANSVERSE  STIFFENERS  C!F  HIGHER  SENDING  STIFFNESS 


In  appendix  A  a  theoretical  solution  for  a  plate  stiffened  by 
stiffeners  of  low  bending  stiffness  is  presented  where  the  buckle 
pattern  was  taken  as  sinusoidal  in  the  longitudinal  direction.  The 
buckle  pattern  of  plates  with  stiffeners  of  higher  bending  stiffness 
is  no  longer  sinusoidal  in  the  longitudinal  direction.  It  is  then 
necessary  to  consider  deflection  functions  which  are  either  symmetric 
or  antisymmetric  about  the  midpoint  of  each  bay  and  are  periodic  over 
an  integral  number  of  bays.  The  critical  shear  stress  of  plates  with 
transverse  stiffeners  of  higher  bending  stiffness  is  analyzed  by  means 
of  the  Lagrangian  multiplier  method. 

Deflection  functions.—  The  correct  buckle  configuration  for  any 
given  plate— stiffener  combination  is  that  which  corresponds  to  the 
lowest  buckling  load.  Several  types  of  configurations  are  investi¬ 
gated.  These  buckling  configurations  are  represented  by  the  following 
two-dimensional  trigonometric  series  (the  coordinates  are  given  in 
fig.  4).  Symmetric  buckling,  periodic  over  each  bay: 


v  =  XI  2 


m=2,4, ...  n=2,4, . 


,  mux  .  njty 
a  sin  — r—  sin  — — 

mn  d  -h 


x:  2 


,  mnx  .  njty 

cos  — sin 
mn  d  b 


(Bla) 


m=0,2, ...  n=l,3. 


Antisymmetric  buckling,  periodic  over  each  bay: 


00 

W  =  5Z 

m=2,4, ... 


x: 

3 3  • •  • 


a 

mn 


sin 


nocx 


sin 


n  Jty 


00 

+  ^ 
m=0,2, . . 


00 

,  1 mac  .  nit  j 

>  cos  — —  sm  — - 

-  mn  a.  Id 

n— 2 , 4, . . . 


(Bib) 
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Symmetric  buckling,,  periodic  over  two  bays: 


I 


w 


m-1 , 3  9  •  •  • 


00 

xi 

n=l,3,... 


a 

ran 


.  rartz  nity 

sin  -  sm  - 

d  b 


00  CD 

+  XI  2Z 

m— 1,  3,  . . .  n — 2,^-, .  •  • 


ran 


cos 


nutx 

d 


sin 


raty 

b 


Antisymmetric  buckling,  periodic  over  two  bays: 


00  00 

w- XT  XI 

m-l, 3, . . .  n— 2,^, • • • 


.  mjtx  .  nity 
amn  sln  “Jf  3in  ~ 


00  00 

+  XI  ^ 

m-l, 3* • • •  -a— 1, 3y • • • 


mjtx  nrty 
b  cos  -  sin  - 

mn  a  b 


Symmetric  buckling,  one  bay;  antisymmetric  buckling,  next 
over  four  bays: 


v-x:  x 


amn 


m=l,3, ...  n=l,3, ... 


m-l 

sin  SS  +  (_i)  2  coa  mrtx 
2d  2d 


XI  X 


mn 


m=l,3,...  n=2,4. 


m-l 

.  mitx  /  t  \  2  nutx 

sin  — - (-1)  cos  - 

2d  2d 


* 


.  (Blc) 


(Bid) 

bay;  periodic 


sin 

b 

ain  22  (B1o) 
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Careful  study  has  shown  that  other  "buckle  patterns  would  require  higher 
buckling  loads  and  that  only  the  five  "buckle  patterns  given  need  be 
considered. 

These  deflection  functions  all  satisfy  term  by  term  the  conditions 
of  simply  supported  edges  at  y  =  0,  b  and  continuity  of  the  plate 
across  the  stiffeners  at  x  -  0,  d,  2d,....  The  condition  that 
stiffener  deflection  equal  plate  deflection  at  the  stiffeners  is 
introduced  by  means  of  Lagrangian  mult'ipliers. 

The  deflection  functions  (Bid)  and  (Ble)  are  found  to  be  the 
governing  ones  for  the  aspect  ratios  investigated;  the  others  lead 
to  unconservative  solutions.  Buckling  criterions  for  the  critical 
shear  stress  are  derived  for  the  deflection  functions  (Bid)  and  (Ble). 

Antisymmetric  buckling,  periodic  oyer  two  bays.—  The  deflection 
of  the  plate  is  given  by  equation  (Bid)  as 


00 


00 


=  21  2 


natx  nfl  y 

a  a  in  — —  sin 
um  d  t) 


n— •  •  • 


00 


oc 


22  2 


mitx  nity 
b,„  cos  — —  sin  — — 
ran  d  t 


m^l, 3 y  • o .  n — 1 y 3, . 


The  deflection  of  the  i*th  stiffener  is  taken  as 


00 


(vs)i  =  2 


Aai  Bln 


n«y 

~b~ 


(B2) 


n— 3^ •  •  • 
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where,  since  the  interval  to  he  considered  includes  two  stiffeners,  i  =  1 
and  2.  The  boundary  conditions  that  stiffener  deflection  equal  plate 
deflection  are 


wCid^)  -  (wJi  =  0  (1=1,2) 


or  upon  substitution. 


+  ^nl  _  0  (n-l,3,...) 

m^l,  3,  * .  • 

bmn  —  Ai2  =  ®  (n=l,3,  •  •  • ) 

IC-  1  y  3  J  *  *  * 


These  equations  show  that  If  A^  is  redefined  as  An 

the  boundary  conditions  become 


bmn+/^  =  0  (n=l,3,...)  (B3) 

m=l,3,... 

These  boundary  conditions  will  be  satisfied  in  the  energy  expression  by 
means  of  Lagrangian  multipliers. 

The  internal  bending  energy  of  the  plate  V,  the  internal  bending 
energy  of  the  stiffeners  Va,  and  the  external  work  of  the  shear 

stresses  T  are  given  by  the  expressions 


Substitution  of  the  deflection  functions  of  the  plate  and  stiffeners 
into  these  energy  integrals  gives 


NACA  TN  No.  I85I 


17 


The  energy  method  requires  that  the  potential  energy  (V  +  Ys  —  T) 

he  minimized  with  respect  to  the  a*s,  h*s,  and  A*s.  Since  the  a’s, 
h*s,  and  A's  are,  however,  hound  hy  equations  (B3),  the  minimization 
is  performed  hy  the  Lagrangian  multiplier  method  hy  minimizing  the 
following  function  F  with  respect  to  the  a*s,  h*s,  and  A’a: 


where  the  7*8  are  the  Lagrangian  multipliers.  When  this  minimization 
is  performed,  the  following  set  of  equations  is  obtained: 


dF 

daw, 

mn 


*  0  = 


^®mn 


16k 


s  h 


00 

Z 


v  mnq. 

VprTjrj 


(m-1,3* • • • ) 

(nS=2j4,  •  e  * ) 


dF 

eft) 


mn 


=  0  = 


V 


mnq 

(^rr° 


l  (B7) 


(m?=l,  3>  •  •  •) 
(n=l,3* • • • ) 


tei 

Dd 


n\ 


+  7 


n 


^-3  3}  •  • . ) 
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When  the  equations  (B7)  are  combined,  the  following  equations  are 
obtained : 


v  ,  2  h2  L  2\2  2EI  4. 

h  m  —  +  nc - n* A, 

mn  l  2  J  Dd  ^ 


'8hsh\2 


m 


2-^m 


jtd  / 

^ .  mr  /  2  2V  2  2' 

\(  2  h2  '  2N2 

q=2j4j . • » 

p-1,3,...  “ftAr 

\  “[2  +  0 

=  0  (B8a) 


0®— 1»3, • • • ) 
(n=*l,3»  •  •  • ) 


Equations  (B8a)  written  in  matrix  form  are 


-  — — 

-  *“• 

— 

Cml  Cml3  Cml5 

V 

**L 

cm31  cm3  Cm35  ••• 

Vi 

=  51 

162A3 . 

Cm51  Cm53  Cm5 

M  •  • 

bm5 

« 

Dd 

1250/^ 

•  •  • 

•  •  • 

• 

• 

• 

• 

(m=l,3,...) 


(B8b) 


where 


S*-f -(¥)“$ 


C  =  vaf 
mn  \ 


mVq2 


q— 2,4, • 


n  —  n 

mnr  mrn 


■■fii 


2  2 

m  q  rn 


=2,4, . . .  (n2  -  q ^(r2  -  q2)(^  ~  +  12^ 
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A  solution  including  all  the  a^'s  and  b^  *s  can  be  obtained 
by  the  first  approximation  of  the  matrix  equation  (B8b) 


Cmlbml  =  2  §  ^ 


(m=l, 3 > • • • )  (B9) 


Substitution  of  V  from  equation  (B9)  into  the  boundary  equation  (B3) 
yields 

0 


The  following  stability  criterion  is  obtained  by  setting  equal  to 
zero  the  coefficient  of 


J.  I 

77  +  —  I  *l  - 0 

Sal 


(BIO) 


m=l,3, ... 


(Bll) 


Similarly,  the  second  approximation  includes  all  the  a  ’s,  b  ,*s 

nm  mi 

and  b  >  *s.  Two  simultaneous  equations  result  from  which  h  and  b 

ml  m3 

can  he  found.  Substitution  of  these  values  into  the  boundary  equation  (B3) 
yields  two  linear  homogeneous  equations  in  A^  and  A^.  If  the  determi¬ 
nant  of  the  coefficients  of  these  two  equations  is  set  equal  to  zero,  the 
following  stability  criterion  is  obtained: 


20 
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Higher  approximations  are  found  in  a  similar  manner.  A  second 
approximation  ■was  found  to  give  satisfactory  results.  For  each  of 
these  approximations,  it  is  necessary  to  find  the  lowest  value 

ET 

of  k  for  each  value  of  — .  The  results  are  given  in  table  I  and 
s  Dd 

in  figure  1. 


Buckling  periodic  over  four  hays.-  The  deflection  of  the  plate  is 
given  by  equation  (Ble)  as 


The  deflection  of  the  i*th  stiffener  is  taken  as 


(B13) 
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vhere  i  =  1,  2,  3,  and  4,  since  the  interval  considered 
stiffeners.  The  "boundary  conditions 


become 


w(id,y)  -  (ws)i  =  0 


m-1 


2 


m=l,3j 


2 


m?=l,  3 3  •  •  • 


V  (-1'  -Aa  -  0 


m-l 

bmn  ^  2  -Aa*° 


m-1 

2 


>  amn  (-1)  2  +  *n2  "  0 


Hk=l,  3j>  •  • 


5 


idfI,  3j  •  •  • 


2 


m-1 

^mn  (“1)  ^  "  ^2  =  ® 


m-1 

®mn  (-1)'5’  *\z‘a 


lb-1 9  3>  • 


> 


hh 1 

bmn  ^  2  +  ^3 


=  0 


m-1,  3  *  •  •  •* 


includes  four 


(i-1,2,3,4) 


(n=lj  3s  *  •  •  ) 


(n=2,4, .  • . ) 


(n— 1, 3^-..) 


(n=2,4, . . « ) 


(n=l,3^...) 


(n-2,4, • . . ) 


22 


NACA  TN  No.  1851 


00  IDe— 1 

y~  =»  (-!>”  -  ^  -  ° 

m— 1, 3  ,  • •  • 

cb  m-1 

y~  v  +\b‘° 

m=l,3.»  ••  • 


(nssl^  3s  •  •  * ) 


(n=2,4j • • . ) 


9 


These  equations  show  that 

Anl  =  _An2  =  ”^3  =  (n-1,3,...) 

Anl  =  Ai2  =  “  A23  =  “  (n=2,4,  . . .) 


If  ^  Is  redefined  as  An,  the  boundary  conditions  became 


m-1 


>  V,  <-«  "A, 


-  A  =  0 


(n— 1, 3 , • • • ) 


m=l,3. 


>  (B14) 


2 


IDr-1 

h  (-1)  2  -  An  =  0 

urn 


(n=2,4, . - . ) 


?b=l-s3s 


These  boundary  conditions  will  be  satisfied  in  the  energy  expression  by 
ire  ana  of  Lagrangian  multipliers. 

The  energy  integrals  are  the  same  as  the  energy  integrals  (B4)> 
except  that  in  the  present  problem  the  upper  limit  of  integration  2d 
is  replaced  by  4d  and  the  upper  limit  of  the  summation  2  is 
replaced  by  4. 


NACA  TN  No.  1851 


23 


The  deflection  functions  of  the  plate  (equation  (Ble))  and  stiffeners 
(equation  (BI3))  are  substituted  into  these  energy  integrals  and  result 
in  the  following  expressions: 


The  minimization  of  (V  +  Vg  —  T)  is  performed  by  the  Lagrangian 

multiplier  method  by  minimizing  the  following  function  F  with  respect 
to  the  a*s,  b*B,  and  A’s. 


where  the  7’ 3  are  the  Lagrangian  multipliers. 


2k 
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When  the  minimization  is  performed,  and  the  resulting  equations 
are  combined,  the  following  set  of  equations  is  obtained: 


CO 


3tt-l 


^mna 


-r  y 


mn  m 


'  b  .  ,  °1_  (-1)  2  =0 

**,*,...  ^(n2-,2)  “  ^ 

(m.— 3^  • « • ) 

(n=lj  3*  •  •  • ) 


\m\rn  rm 


a“1  (  a2  -  n2  ) 


nq 


q=l,3,.. 


m-1 

+  2  E  „V  (-1)  2  =  0 


> 


M 


(m-1, 3, .  •  • ) 

(n=2,h, . . . ) 


where 


J  mV 

A  =  2  — —  +  n 

V  4d! 


8k  b  111-1 

%  -  “(-D 


(Bl6a) 
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Equations  (Bl6a)  in  matrix  form  are 


Substitution  of  a^  and  b^  from  equations  (B17)  into  the  boundary 
equations  (Bl4)  yields 


26 
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4  2 

AmlAm2  —  9  rm 


VlAn 2  “ 


>  (BIB) 


If  the  determinant  of  the  coefficients  of  the  linear  homogeneous 
equations  (Bl8)  is  set  equal  to  zero,  the  following  stability  criterion 
is  obtained: 


(B19) 


Similarly,  from  the  second  approximation,  including  all  the  a^, 

b  and  a  _  terms,  the  following  stability  criterion  is  obtained: 
m2*  m3 
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TABLE  I.  -  SHEAR-STRESS  COEFFICIENTS  FOR  STIFFENED  PLATES 
WITH  PANEL  ASPECT  RATIOS  OF  ONE,  TWO,  AND  FIVE 


PlateB  with  stiffeners  of  low 
bending  stiffness6 

Plates  with  stiffeners  of  higher 
bending  stiffness 

El 

Dd 

kg 

Aspect 

ratio 

Antisymmetric 
buckling 
periodic  over 
two  bays 

Buckling 

periodic 

over 

four  bays 

11 

Dd 

kB 

El 

Dd 

k 

s 

0 

2 

5 

20 

50 

100 

200 

5-3^ 

10.34 

16.07 

37.14 

68.99 

112.2 

184.6 

1 

0 

2.91 

7.78 

22.29 

00 

5-53 

7.85 

9.82 

UL.78 

13.86 

0 

7.09 

19.03 

00 

6.08 

10.0 

10.5 

10.86 

2 

.  ....  -  .... 

0 

3-35 

14.50 

22.99 

33.11 

45-77 

61.97 

82.92 

112.3 

605 

00 

9.65 

12.0 

16.0 

18.0 

20.0 

22.0 

24.0 

26.0 

28.0 

35.0 

37.05 

0 

5.475 

11.93 

26.37 

36.29 

68.92 

145.4 

625 

00 

5.54 

15.0 

20.0 

23.O 

24.5 

26.0 

27.0 

28.0 

28.2 

5 

0 

18.02 

90.99 

176.8 

444.7 

704.4 

42.5 

70 

90 

100 

120 

140 

0 

49.19 

112.8 

220 

00 

agH 

independent  of  aspect  ratio. 

All  the  aml,  bj^,  aj^,  &nd  bj^ij. 
CA11  the  am3,  b^,  a^,  and  b^g 


coefficients  used, 
coefficients  used. 
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TABLE  II .  -  EXPERIMENTAL  BUCKLING  DATA  OF  SHEAR  WEBS 
WITH  UPRIGHTS  NOT  CONNECTED  TO  THE  FLANGES 


Specimen 

d 

b 

t 

Uprights 
(nominal  size) 
(in. ) 

El 

T 

ks 

(a) 

(in.) 

(in.) 

(in.) 

Dd 

(ksi) 

2-D-O 

5-0 

23.5 

0.0397 

l/2  X  l/2  x  l/l6 

221 

2.66 

101 

3-D-O 

5.0 

23.5 

.0394 

3/4  x  3/4  x  1/16 

680 

3.08 

116.5 

4-D-O 

5-0 

23.5 

.0405 

3/4  x  3/4  x  3/32 

946 

3.295 

U7.5 

5-D-O 

10.0 

23.5 

.0404 

1/2  X  1/2  X  l/l6 

98.3 

1.21 

43.3 

6-D-O 

10.0 

23.5 

.04o8 

3/4  x  3/4  x  l/l6 

306 

1.54 

54.2 

7-D-O 

10.0 

23.5 

.0410 

3/4  x  3/4  x  3/32 

456 

1.47 

51.3 

8-S-O 

5.0 

23.5 

.0394 

1/2  x  1/2  x  0.064 

95-8 

2.895 

109 

9-s-o 

5.0 

23.5 

.0399 

3/4  x  3/4  x  3/32 

456 

3.01 

111 

10-S-0 

10.0 

23.5 

.0410 

1/2  x  1/2  x  1/16 

41.4 

.82 

28.6 

ll-S-0 

10.0 

23.5 

.0398 

3/4  x  3/4  x  1/16 

151.5 

1.357 

50.1 

12-S-O 

10.0 

23.5 

.0405 

3/4  x  3/4  x  3/32 

217 

1.41 

50.3 

lSj  stiffeners  on  one  side  of  plate. 

By  stiffeners  on  loth  sides  of  plate. 
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TABLE  III .  -  EXPERIMENTAL  BUCKLING  DATA  OF  THICK  WEB  BEAMS 
WITH  UPRIGHTS  CONNECTED  TO  THE  FLANGES 


Uprights 


Specimen 

(a) 

L 

(in.) 

d. 

(in.) 

Y-12-7S 

9.88 

7.00 

V— 12— 8s 

9.88 

7.00 

V-12-9D 

9.13 

7.00 

Y— 12— 10S 

' 

9.88 

7.00 

V— 12— 11D 

9.13 

7.00 

V— 12— 12S 

9.88 

7.00 

V— 12— 13D 

9.13 

7.00 

V— 12— 14S 

9.88 

7.00 

Y— 12— 15D 

9.13 

7.00 

(nominal  size)  M  (]fal) 


0.1005 

.1044 

.1025 

.1043 

.1025 

.0987 

.1000 

.1007 

.1057 


(in.) 


li  X  li  X  i 

8  8  8 


91.0  15.5  15.4 


3/4  x  3/4  x  1/8  25.8  15.4  14.15 

5/8  x  5/8  x  1/8  4o.4  16.8  13.65 

5/8  x  5/8  x  1/8  14.5  16.3  15,0 

5/8  x  5/8  x  3/32  30.3  17.2  14.0 

1/2  x  1/2  x  1/16  4.1  12.3  12.7 

1/2  x  1/2  x  1/16  11.3  13.1  11.15 

5/8  x  5/8  x  3/32  11.2  13.2  13.1 

5/8  x  5/8  x  1/16  18.8  15.7  12.O 


S,  stiffeners  on  one  side  of  plate. 
D,  stiffeners  on  "both  sides  of  plate 
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Figure  1.-  Continued. 
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Figure  3.-  Thick  web  beam  (with  uprights  connected  to  flanges)  after  failure. 


